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FOREWORD 
The work described in this report was performed between January 1971 and January 1972 
as a part of Project 2567 in support of munition development programs such as the Close Air 
Support Weapon and the Modular Weapon Series at the Air Force Armament Laboratory, 
Eglin Air Force Base, Florida. 
This technical report has been reviewed and ts approved. 
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ABSTRACT 


A computer program using an iterative technique was developed for determining the radius 
of action of F-4E and A-7D aircraft with any configuration of external stores. The program 
uses performance data from the aircraft flight manuals to calculate fuel and distance required 
to achieve military power climbs and optimum cruise fuel consumption. Required inputs to 
the program include the initial amount of fuel onboard and aircraft gross weight, the cruise 
and loiter altitudes, and the outbound and returning drag indicies due to aerodynamic drag 
of the external stores. Optional inputs that provide increased accuracy include fuel for 
engine start and taxi, and fuel and distance for takeoff and descent. The program calculates 
sequential increments of the outbound cruise portion of the mission profile until the fuel 
reserve desired at the end of mission (return to home base) is obtained; thus the radius of action 
is determined. 


Distribution limited to U. S. Government agencies only; 
this report documents test and evaluation; distribution 
limitation applied January 1973. Other requests for this 
document must be referred to the Air Force Armament 
Laboratory (OLGC), Eglin Air Force Base, Florida 32542. 
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SECTION | 
INTRODUCTION 


Carriage of external stores (weapons, fuel tanks, instrumentation packages, and suspension 
equipment) can result in large degradations to the range and speed performance of tactical 
fighter aircraft such as the F-4E and A-7D. Some weapon configurations exhibit aerodynamic 
drag characteristics that are as great as those of the aircraft itself; consequently, the aircraft 
range performance may be reduced as much as 50 percent. Therefore degradation plays a 
significant role in overall mission effectiveness, and the store developer must consider aerodynamic 
drag as an important parameter in the many trade-offs leading to the final design. This report 
describes a computer program developed for evaluating range performance effectiveness. The 
program can readily be used on a parametric basis and will provide a quantitative method 
of assessing the effect of aerodynamic drag on the aircraft combat radius of action. 


The validity of the program results depends primartly on the accuracy of the store drag 
input data. Both the F-4E and A-7D arrcraft utilize a drag index system based on the 
increment of aerodynamic drag due to external stores. Each external store authorized for 
carriage has a drag index assigned which included interference effects resulting from the 
additional drag of stores in close proximity to one another. The total drag index of a 
particular configuration is obtained by summing the individual drag indicies of the external 
stores as described in the appropriate aircraft flight manuals (e.g., References 1 and 2). 
Precise drag indicies cannot be determined without extensive wind tunnel tests and/or flight 
testing; however, good approximations for an untested store can be obtained from wind tunnel 
data of similar stores. References 3 and 4 describe techniques for the latter procedure, and 
methods for reducing wind tunnel data to obtain more accurate drag indicies may be found 
in Reference 5. A typical mission profile ts shown in Figure 1. Performance data from the 
F-4E and A-7D aircraft flight manuals (References 1 and 2) are used to calculate the fuel 
and distance required for military power climbs and optimum cruise fuel consumption rates. 
Inputs to the program are the initial fuel weight and aircraft gross weight, the cruise and 
loiter altitudes, ar-d the outbound and returning drag indicies of the external stores. The amount 
of fuel required for engine start, tax:, takeoff, and descent, as well as the distance required 
tor takeoff and for descent, are optional inputs that can provide increased accuracy. The 
program caiculates sequential increments of the outbound cruise portion of the «mission 
profile unti! the fuel reserve desired after returning to base is obtained; thus determining the 
radius of action. The program can be used for any mission profile provided that the cruise 
portion is conducted at the optimum condition for the altitude selected and that a military 
power climb to altitude is assumed. 
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SECTION I! 


DEVELOPMENT OF THE ITERATIVE TECHNIQUE 


The fuel and distance required for the climb and cruise portions of the mission are highly 
dependent upon the outbound and returning altitudes and drag indicies; however, the fuel 
and distance for engine start, taxi, and takeoff are independent of these parameters. Also, 
a good approximation for descent from any altitude for any drag index is that the fuel 
expended and the distance traveled are both zero. Therefore, for parametric studies, engine 
start, taxi, takeoff, and descent can be considered as constants. Guide lines for obtaining 
accurate estimates for these portions of the mission are contained in the aircraft flight manuals. 
The program is designed to accept the following estimates as input data: 


a. 


b. 


Cc. 


d. 


e. 


Fuel and distance traveled away from home base for start and taxi. 
Fuel and distance traveled away from home base for takeoff. 

Fuel a..d distance traveled away from home base for outbound descent. 
Fuel and distance traveled away from target for descent returning. 


Fuel and distance traveled away from home base during loiter in target area. 


Additional required program inputs are as follows: 


a. 


b. 


g. 


h. 


Aircraft type. 

Altitude outbound in the target area, and returning. 

Fuel at engine start. 

Payload (any weight which might be dropped, ejected, or jettisoned during mission). 
Basic aircraft weight (gross weight minus fuel and payload). 

Fuel reserve at end of mission. 

Payioad not dropped, ejected, or jettisoned during mission. 


Drag index outbound and returning. 


The drag indicies outbound and returning from the target area should be determined by 
the best available means. 


For each set of input data, radi of action are determined for loiter fuels of 9600 pounds 
down to 0 in increments of 12090 pounds. Obviously, the amount of fuel available for loiter 
can be used to approximate a permissible loiter time. 


If the altitude in the target area is greater than zero but less than 5000 feet, th > nregram 
automatically changes to either zero or 5000 feet based upon the magnitude of the input 
vatue. This conversion is necessary because the subroutines are unabie to determine accurate 
values for this range of altitudes. 


Further, if the aircraft gross weight at engine start exceeds that allowab'e for the aircraft, 
fuel is downloaded automatically in the program to the point where the gross weight is within 
allowable limits. 


The weight and distance traveled away from home base ts determined by subtracting all 
fuel used and adding all distances traveled. Subroutines are then called which provide the 
fuel and distance required to climb to the outbound altitude; at that time new weight and 
distance values are determined and retained. If this is followed by a cruise, it will be a 
1000-pound-fuel cruise, and the gross weight at the beginning of the crutse will be used. 
The estimated subsequent fuel requirements determine whether a 1000-pound-fuel cruise 
is undertaken. Another subroutine provides the fuel consumption rate. After the outbound 
climb or the 1000-pound-fuel cruise, as the case may be, the outbound descent 1s made, a 
9600-pound-fuel loiter is completed; all, none, or some portion of the payload ts dropped; 
and the return trip home is begun after the distance from home base to the target 1s noted. 


Each of the climb subroutines is called twice to determine the fuel and dstance required 
to proceed from target altitude to the return altitude. The amount required to climb from sea 
level to the target altitude is subtracted from the amount required to climb from sea level 
to the return altitude. 


After the return climb 1s completed, the cruise subroutine 1s called to determine a fuel 
consumption rate for the beginning of the return cruise. With the fuel data for the return 
descent and for landing and taxi, the desired gross weight at the end of the return cruise and 
the fuel consumption rate at this point are Known. The same distance data, the return climb 
distance, and the distance from home base to the target, are used to calculate the return 
cruise distance. The fuel consumption rates are averaged, and the fuel required for the return 
cruise is computed. 


The trial mission 1s completed and the extra fuel determined. If the extra fuel ts less 
than minus one pound, the mission can not be completed for the specified loiter fuel, and 
another simulated mission 1s attempted with 1200 pounds less fuel for Joiter. If the extra 
fuel amounts to one pound or more, either one-half, one-third, or 1000 pounds of it are 
used for a cruise or an additional outbound cruise to the target area. The appropriate 
quantity to be used is determined by the relation of the outbound and returning drag indices 
combined with the amount of extra fuel. The iteration is repeated using either the weight 
and distance after the outbound climb or the weight and distance after the last outbound 
cruise, if one has occurred, until the mission is completed with the desired quantity of extra 
fuel (zero or minus one pound). Thus, the radius of action has been determined. 


Subsequent radii are determined for loiter fuels incremented by minus 1200 pounds until 
the fuel for loiter has been reduced to zero pounds. Instead of simulating the new mission 
from the very beginning, the new iteration is begun after the outbound climb or the last 
1000-pound-fuel cruise, if one has ¢ :curred. 


Appendixes |, 11, and #1} contain the radius-of-action program listing, input list, and 
limitations, respectively. 
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SECTION III 
DEVELOPMENT OF THE SUBROUTINES 


This section discusses the development of the subroutines for calculating the fuel and distance 
required to climb and the fuel consumption rate, given the aircraft type, altitude, gross weight, 
and drag index. 


The F-4E and A-7D aircraft climb and cruise parameters are displayed in figures contained 
in the aircraft flight manuals (References 1 and 2). These figures define the fuel consumotion 
rate for optimum cruise, fuel required te climb, and distance required to climb. These figures 
are set up in a parametric format and are a function of aircraft gross weight, altitude, drag 
index, and a parametric variable (designated as T in this report). Each figure is divided into 
two charts with the first chart being a family of curves that are a function of gross weight 
and altitude which define the dependent variable T. The second chart has a family of curves 
with independent variables of T and aircraft drag index and with a dependent variable of either 
fuel for climb, distance for climb, or fuel consumption rate for optimum cruise depending upon 
which figure is being used. Although no equations define the families of curves, a method 
was developed to define the families for computer applications. The first chart contains 
discrete curves which are a function of gross weight, for altitudes of 0, 5000, 1000, .. ., and 
40,000 feet. Similarly, the second chart contains discrete curves for drag indicies of zero to 
the maximum in specified increments that are a function of T. The problem then was to find 
equations of lines for any condition other than those discretely defined. The method used 
was a Newtonian interpretaticn technique (Reference 6) that defines an nth order polynomiat 
for n discrete points. The development of the Newtonian equation by this methoa forces 
the polynominal to satisfy the n discrete points by definition. The equation then can be used 
to calculate a value between ai-y of the discrete points. 


it was hypothesized that if a family of straight parametric lines was only defined at discrete 
values, the slope and intercept of any line lying between two defined lines could be determined by 
interpolation of the defined values. Although the family of curves did not consist of straight 
lines, the values were approximated by decomposing the curves into a number of straight lines. 
In the first chart, the curves that are a function of gross weight and altitude are drawn as straight 
lines between certain gross weights. Similarly in the second chart, the curves that are a function 
of T and drag index are drawn as straight Itnes between certain values of T. Now using the 
Newtonian method, two equations are defined to give the slope and intercept for any altitude 
line in the first chart. Similarly, two equations are defined that provide the slope and intercept 
for any drag index line in the second chart. Thus, when the altitude ts defined, a linear equation 
which has gross weight as the independent variable and T as the dependent vartable, 1s con- 
structed with a calculated slope and intercept. When a drag indew is defined, another linear 
equation, which has T as the independent variable and either fuel for climb, distance for 
climb, or fuel for cruise as the dependent variable, is constructed with a calculated slope and 
intercept. 


Three computer subroutines were constructed to incorporate the method just described. 
One subroutine will construct the equations for determining the fuel for climb, a second 
subroutine will give the equation for distance to climb, and a third subroutine is concerned 
with the fuel consumption rate for optimum cruise. 
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Appendixes IV and V contam the main programs and two additional subroutines for 
utilizing the three radius-of-action subroutines to obtain tabular data for the following 
independent variables: 


a. 


F-4E (see Appendix V1) 

(1) Altitude (5000 to 40000 feet in increments of 2500 feet). 

(2) Weight (30000 to 58000 pounds in increments of 1000 pounds). 
(3) Drag index (0 to 140 in increments of 10). 

A-7D 

(1) Altitude (5000 to 40000 feet in inc. 2nents of 2500 feet). 

(2) Weight (20000 to 42000 pounds in increments of 1000 pounds). 


(3) Drag index (0 to 300 in increments of 20). 
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SECTION iV 


APPLICATION OF THE TECHNIQUE 


As an example, assume that the radius of action for an F-4E aircraft is to be determined 
for a high-low-high mission at altitudes of 20000 feet outbound, 5000 feet in the target area 
(loiter), and 30,000 feet returning. Further, assume some realistic values for the mission 
segments (Figure 1) as follows: 


Fuel, tb Distance, N. M. 
Start and Taxi 440 0.5 
Takeoff 925 1.5 
Descent going to target 300 10.0 
Descent returning from target 200 15.9 
Landing and taxi 100 2.0 
Loiter N/A 0.0 


Also, assume that the payload ts ssx 800-pound weapons; all of which are to be dropped 
from the triple ejector rack on the inboard armament stations. A typical fuel reserve at the 
end of a mission is 2000 pounds. 


The aircraft incremental drag as a result of weapon carriage (weapon drag index) 1s largely 
dependent upon weapon shape and diameter and can conceivably vary from a total weapon 
drag index of 30 to 60. Figure 2 shows the radius of action as a function of loiter fuel for 
two weapon drag indicies and three fuel loadings. Note that when fuel was carried on the 
centerline and outboard stations with the large payload, 1198 pounds had to be downloaded 
sO as not to exceed the maximum gross weight of 58000 pounds. This makes the increase 

in the radius of action based on outboard and centerline fuel somewhat less than if the 
payload was less. 
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Figure 2. Radius of Action Versus Loiter Fuel for Two Weapon Drag Indices and Three 


Fuel Loadings 


APPENDIX | 


RADIUS OF ACTION PROGRAM LISTING 
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PROGRAM P3700C INPUT, OUTPUT, TAPES=INPUT) 

DIMENSION OACR(100) ,O0TT(100) 

INTEGER PAW, F,FAS, FECL) FFCR FFOG yFFOR,FFLyFFLAT yFFSAT yFFTOsFTBRyPy 
1PE,PRyW,WACR( 100) sNAEOCR yXFFCRyWALBCRy XFyFFCLL yFFCL25C 

802 FORMAT(A10,3F5.0,54X , 11) 

COMMON/CACFT/F4E,A70 

802 FORMAT (BA10) 

B03 FORMAT (314 ,4X 4TH, 2F 3004 2X 9215916 6F 400) 

B01 FORMAT 31X%y20T6y4X) y3Xy IT3y SXeEC F601 94) 9/9 1K, SUITE 4X) 9 3Ky 13 
195%, 7(T6 4X) 9 /) 

902  FCRMAT(1H1,42Xy GHFUEL 7) 32Xy7HPAYLCADy 2X y BHRET URNEDy 5X) 4HORAG » 5X95 
LHSTART,2013Xy THOESCENT) 9 4X y 7HLANDING y/y 2X y SHBASIC ys 5X ySHTOTAL 9 1 4X 98 
Z2HRETURNED, 1X» 1OHPROGRAMMED, 3X4 SHINDEXy 3X, B8HAND TAXI y2X,8HTAKE-OFF, 
33X,5HGCING, 5X, 6HLOL TER, 2X, SHRETURNING, 2X5 8HAND TAXI,4X,6HRADIUS,/y 
42X,4HACET, 6X5 SHFUEL » 5X» 7HPAYL OAD y 3X y BHEXPENOED y 3X4 OHACTUAL y 3X y SHRE 
STUPNINGs 4X, 2HNMy5 (BX p2HNM) 9 5Xy9HCF ACTION)/91X 9 2(6HWT (LB) 4X) 1X93 
6 (BHNT (LB) 9 4X) 9 1 XySHGCING 6 (4X, EHWT (LED) 9 7%y2HNMy/) 

903  FORMAT(SiX,*CHECK THE INPUT DATA FOR THIS CASEe*,/) 

904 FORMAT(///,24X,*THIS PROGRAM DOES NOT PRESENTLY HAVE REQUIRED INFO 
4PMATICN FOR THE *,A10,* AIRCRAFT,.®) 

905  FORMAT(///,47X,*ALL ALTITUDES MUST BE SEA LEVEL OR ABOVE.®) 

906 FORMAT(///,31X,*AN ALTITUDE OF *, I5 ,* FEET EXCEFOS THE TABLES F 
10R THE *,A10,* ATRCRAFT.*) 

907 FORMAT(29X,*A GROSS WEIGHT OF *,15,* POUNOS EXCEEDS THE TABLES FCR 
4 THE *,A10,% AIRCRAFT.*,/) 

908  FORMAT(34X,*A ORAG INDEX OF *, I4 4* EXCEEOS THE TABLES FOR THE * 
1,A810,*% AIRCRAFT. *%,/) 

909 FORMAT (1H1,55X,*RADIUS OF ACTION FOR THE*%s//,59X yA10y *AIRCRAFI% 4 // 
4,49X,*ALTITUDE GOING TO THE TARGET="%, I5, * FEET*,//,46X,*ALTITUD 
2E RETURNIAG FROM THE TARGET-*,15, * FEET%,»//,50X,*ALTITUDE AT THE 
3 TARGET-*,15, * FEET®) 

910 FORMAT(1X,*A GROSS WEIGHT OF *,I16,* POUNDS AND A ORAG INOEX OF *, 
413,* AND AN ALTITUDE OF *,15,* FEET EXCEEO THE TABLES FOR THE *, 
2A10»* AIRCRAFT. %,/) 

G11 FORMATO 31X92 0T6y4X) 9 2Xy Fe Dg 5X96 FG0t 4X) oy 1X pS L6y 4X) 42K y FHL 
195Xy6(16,4X)) 

912 FORMAT (31X,20T6y4X) y3Xy T3y SX 60 F601 ,4X) Ay 1X yp 5(16,4X) y3Ky 13 
195Xy7(16 4X) 9 /y1H1) 

913 FORMAT(1X, G2X,4HFUEL»/y 32Xy THPAYLOADs 2X5 BHRETURNED, 5X y 4HDRAGy 5X95 
LHSTART,2(13Xy THDESCENT) 4 4X y PHLANDING 9/4 2X ySHBASICy 5X SHTOTAL y 14X98 
2HRETURNED, 1X 1OHFROGRAMMED y 3X SHINDEX, 3Xy8HAND TAXI, 2X, 8HTAKE=OFF, 
33Xy SHGOINGy 5Xy GHLOI TER, 2X, SHRE TURNING, 2Xy BHAND TAXI 4X ,GHRANTUS, 7, 
4 2X, GHACFT ,6X y GHFUEL y 5X4 7HPAYL OAD y 3X, BHEXPENDED » 3X y GHACTUAL » 3X5 QHRE 
S TURNING 9 4X y2HNMy5 (8Xy2HNM) 5X9 SHOF ACTION y/ 1X5 2(6HWT (LB) 4X) 51X53 
6 (GHWT (LB) 9 GX) 9 1X9 SHGOING 96 (4X, EHWT (LB) ) » 7Xy2HNMy/) 

914  FORMAT(1H1,*ALL WEIGHTS ARE IN POUNOS. ALL OISTANCES ARE IN NAUTI 
1CAL MILES, ALL FUEL CONSUMPTICN RATES EXCEPT THE AVERAGES ARE NAU 
2TICAL MILES PER*y/, 
31X,*POUNN. AVERAGES ARE POUNDS PER NAUTICAL MILEe*)//////y 
G1iX,*WEIGHT AT THE ENO OF RETURN CRUISE=*,16,*% FUEL RATE=*sF 504, // 


5/) 

g15 FORMAT(1X,*FOR CLIMB GOING TO THE TARGET WEIGHT=*,16,* FUEL=* 
teIG,* OQISTANCE=*%,F7.4y/) 

916 FORMAT(1X,*FOR CRUISE NUMBER*,I3,* GOING TO THE TARGET WEIGHT= 


L*sl6y/, 


11X,*FUEL RATE=*,F5.4,*% TOTAL DISTANCE TRAVELED AFTER THIS CRUISE= 
2*,FB eu y/) 

917 FORMAT(1X,*FOR CLIMB RETURNING FROM THE TARGET WEIGHT=*%, 1695/5 
16X¥FROM 0 FEET TO RETURN ALTITUDE FUEL=*I4&* OISTANCE=",F7.4/ 
1LOX*FRCM 0 FEET TO TARGET ALTITUDE FUELS=*IG* OISTANCE=%,F7.4/ 
L6X*FROM TARGET TO RETURN ALTITUDE FUELH*IG*® OISTANCE=*,F7.4/ 
41) 

918 FORMAT (1X,*WEIGHT AT BEGINNING OF RETURN CRUISE=*%,16,% FUEL RATE 
4 AT THIS POINT=%,F5.4,% AVERAGE RATE=",F8.4,/; 
24X,*DISTANCE FOR RETURN CRUISES", F8.4,*% FUEL AT ENO OF MISSTONS*, 
3I5,* EXTRA FUEL AT ENO OF MISSION=*%,15,////) 

919 FORMAT(1X,*FOR CRUISE NUMBER*,I3,* GCING TO THE TARGET WEIGHT= 
1*,T6e/,y 
LIX,* FUEL RATE=*,F5.4,% EXTRA FUEL FCR CRUISE GOING=*y5I5y5* TOTAL 
20ISTANCE TRAVELED AFTER THIS CRUISE=*%,F8.4,/) 

10 ASSIGN 212 TO IFEOF 
! CALLECF(TEOF) 
FuE=10HFGE 
A7C=10KA70 
READN(5,801) ACFT,ALTG,ALTRAL TAT ,C 
IFCALTAT LES C500 ee ANCAALTAT.GTo 0.) ALTAT=0. 
TF CAL TAT LT eo S000. eANDLALTAT.GT.2500.) ALTAT25000. 
ILTG=ALTG 
ILTR=ALTR 
: ILTAS=ALTAT 
| PRINT 90¢, ACFT,ILTG,ILTR,ILTA 
TF CALTGeLTeDe cOReALTReLToDe ORsALTATs6LT 602) GO TO 202 
| ILT=ALTG 
IF (CALTG.LT.5000.) GO TO 203 
ILT=ALIR 

| TF CALTR.LT.50006) GO TO 203 

| ILT=ALTG 
IFCALTR.GT.ILT 9 ILT=ALTR 
IF (ACFTCEQ.FUE. AND ILTeGTe40000.) GO TO 203 
IF CACFET.EQeA7TD AND (ILTGT»4O090-ORILT oLT.5000)) GO TO 203 
IF CACFT LEQ eFYESOReACFTEQ.A70) GO TO 4&0 

20 PRINT 904,ACFT 

2t ASSIGN 211 TO IE0F 
00 30 I=1,10000 

306 READ(5,802) A 
GO TO 212 

40 ASSIGN 211 TO IEOF 

50 IF(CeNE1) PRINT 902 
00 200 IP=1,2 
READ(5,803) FFSAT,FFTO,FFOG, FFOR,FFLATyFIBRyPRUIGeDIR, 

LFAS,P,BAW,OFSAT,OFTO,OF0OG,OFLyOFCR,OFLAT 


Bad 


ICR=CIP 
IO0G=01G 
FFL=9€00 
W=RAW+FAS+P 
IF CATE TCEQeF HES ANDeHeGT258000) FAS=58000-BAW-P 
TF CACFTSEGsA7TN AND eWeGTe 42000) FAS=42000~BAH-P 
IFCFFOAT L120) GO TO 201 
IF(FriC «LT.0) GO TO 201 

| IFCFFOG .LT.0) GO TO 201 

H 

i 

i 

i 
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i IFCFFOR .tT.0) GO TO 201 
‘ IFCFFLAT .LT.0) GO TO 201 
IFCFTAR .LT.0) GO TO 201 
IF (FAS eLE.0) GC TO 201 
IF (PR LT.0)? GC TO 201 
IF CP eLT20) GO TO 201 
TF (BAW LE.0) GO TO 201 
IF (FAS. LEC FFSAT+CFFI CES FOG +FFOR+FFLATSFIBR)) GO TO 201 
IF(PR.GT.F) GO TO 201 
, PE=P-PR 
J=2 
10=N IG 
IFCOIR,LT.IN) IO=NIR 
: TF(IOeLT.06) GO TO 205 
10=01G 
IF(DIR.GT.IN) IO=OIR 
IF CACFTSEQeF UE ANDO. IDoGT0140 ) GO TO 205 
IF CAGFT.EQ.A7D.AND.ID2GT 2300 ? GO TO 205 
IFCOIR,GT.OIG) Js3 
XJ=.! 
Jamss*1000 
WAEOCR=BAN+PREFTBREFFLAT+F FOR 
IF CACFTEQsF4E AND WAEOCR.LT. 30000) GO TO 210 
IF CACFTSEQ.A7N. ANDeWAEOCR»LT.20000) GO TO 210 
WsBAW+FAS4P 
W=W-FFSAT-FFTO 
ILT=ALTR 
I0=D0IR 
CALL CLOCACFT,ALTR,HAEOCR,DIR,FRE), RETURNS (207) 
ILT=ALIG 
IO=0 IG 
CALL CROCACFT,ALT&,WAEOCR,DIR,FRE) 
TF(C6EQ.1) PRINT 914, WAEOCR, FRE 
O=OF SAT+0FTO 
CALL CLF CACFT,ALTGyW,OIGyFFCL) 
CALL CLOCACFT,ALTG,WyDIG,OFCL }y RETURNS(207) 
IFFCL=FFCL 
IF(CeEQ.1) PRINT 915, W,FFCL ,OFCL 
W=W-FFCL 
O=O0+0FCL 
I=2 
IMO=I-1 
WAGR(1)=W 
DACR(1)=0 
IAs! 
i NXFFCR=3001 
NOI=9 
00 199 IFFL=1,NOI 
IF (CF AS@FFSAT-FFTO-2* IF FCL<FFOG-FFL@FFOR=FFLAT“FTBR).LTe JBM ) GOT 
10 101 
IF(NXFFCR.LT.1000) GO TO 1014 
I=IA 
IMO=I-1 
106 CALL CROCACFT,ALTG,WACR(IMO),0IG,0FCR) 
WACR(I) =WACRCIMO)=1000 


' DACR(T)=D0ACR( IMO) 40FCR*10006 
: IF(CeEC.1) PRINT 916, IMO,WACR(IMO) ,DFCR,DACR(T) 
WALACR=HACR(T) 
DALACR=DACRIT) 
T=I+t 
TAzI 
IMO=T-1 
404 W=WACR(TVO)-FFOG-FFL<PE 
OTTCIMC) =DACRCIMO) OF DG 40FL 
CALL CLF CACFT,ALTRy»W,OIR,FFCL2) 
CALL CLE CACFT,ALTAT,W,DIR,FFCL1) 
CALL CLOCACFT,ALTR,W,OIR,DFCL2), RETURNS (207) 
CALL CLOCACFT,ALTAT »W,OIR,OFCL1) »RETURNS (207) 
FFCL =FFCL2-FFCLI 
OFCL=CFCL2-OFCLI 
TF(C.EQe1) POINT G17, WyFFCL2,0FCL2,FFCL1,0FCL 1 ,FFCL,OFCL 
W=oW-FFCL 
DFCR=DTT (IMO) -DFCL-OFOR=0FLAT 
CALL CFOCACFT, ALTR yW,OIR,FPRB) 
| FRECQ=2./(FREFFRA? 
‘ FFCR=FFFCR*DFCRE.S 
F=WeBAR-PR-FFCR@FFOR=FFLAT 
XFrFeFTRR 
TF(C.EQe1) PRINT 9185 WeFRBgFRFECR,OFCRy, Fy XF 
KFFCR=XF/X5N4.5 
IF (XFFCR.GE.1000) GO TO 100 
IF (XFFCRoLTe-1) FFL=FFL=1200 
IFCFFLeLT.0O)2 GO TO 200 
IF(XFFCR.LT.-1) GO TO 101 
IF(XFFCReLT.1) GO TO 198 
NXFFCREXFFCR+1200/J 
| 402 CALL CROCACET,ALTG, WACRCIMO)»DIG,FR) 


os 


WACR(T) =WACRC IMO) -XFFCR 
DACR(T)=DACRCI IMO) CF REXFFCR 
TF(CeECe1) PRINT 919, IMO,WACROCIMO) ,FR,XFFCR,OACR(T) 
Ts=I+1 
IMC=T-1 
W=WACR (IMO) <FFOG-FFL SPE 
OTT( IMC) =CACR(IMC) +0FOG+OFL 
CALL CLF CACFT,ALTRyW,DIR,FFCL2) 
CALL CLF CACFT,ALTAT,W,OIR,FFCL1) 
CALL CLOCACFT ,ALTRyW,OIR,DFCL2), RETURNS (267) 
CALL CLOCACFT,ALTAT WyOIR,DFCL1), RETURNS (207) 
FFCL=FFCL2-FFCL1L 
OFCL=CFCL2-OFCL1I . 
IF(Ce&Q.1) PRINT 917, WyFFCL2,0FCL2,FFCL1,0FCLAi,FFCL,OFCL 
WHW-FFCL 
DFCR=CTT CIEMO) -OFCL=OFOR=DFLAT 
CALL CROCACFT,ALTRyW,OIR,FRB) 
FRFCR=2,/{(FRE+F RB) 
FFCREFRFCRYDFCR4+.S 
F aW-BAK-PR-FFCR-FFOR<FFLAT 
XFaF-FTBR 
TE (C.ECe1) PRINT $18, WyFRE,FRFECR,DFCR»F yXF 
XFFCR=XF/XJS+05 


: a a RE 
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198 


199 
200 


212 


IF(XFFCR.GT.20) GO TO 102 

IOT=OTT(IHO) +25 

TF(C2EQ.1) PRINT 913 

IF(CoNEe1) PRINT 901, PR,»FTBR,IOR,DFSAT,OFTO,OFOG,OFL,OFOR,DFLAT 4B 


LAW FAS, P, FEF TOG sFFSAT,FFIC,FFOG FFL »FFOR,FFLAT,IOY 


IF (CeEQe1) PRINT 912, PReFTBR,ICRyOFSAT,OFTO,OFOG,DFLyDFOR,DFLAT ,B 


LAW,FAS,P PE, Fy ING sFFSAT,FFIO,FFOG,FFLyFFOR,FFLAT,IDT 


FFL=FFL=1200 

IF(FFL.LT.0) GO TO 200 

I=IA 

IMO=T*4 

CONTINUE 

CONTINUE 

GO TO 50 

CONTINUE 

PRINT 914, PRyFTBR,DIR,OFSAT,DFTO,DFDG,DFL»DFOR,DFLAT » BAN, FAS» Py PE 


LyoFsDIGyFFSAT,FFTO,FFOG,FFL,FFOR,FFLAT 


PRINT 903 

GO TO 200 

PRINT 905 

GO TO 21 

PRINT 906, ILT, ACFT 
GO TO 24 

H=4HFFCL 

PRINT 907, ALT yW Hy ACFT 
GO TO 200 

PRINT 908, IO, ACFT 

GO TO 200 

PRINT O10, W,IO,ILT,ACFT 
GO Te 200 

PRINT 907, WsACFT 

GO TO 200 

GO TO 10 

CALL EXTT 

ENO 


SUBRGUTINE CLOCACFT,ALT, W,DI,0IS), RETURNS (A) 
DIMENSION DY5(10) ,0Y6(10) ,DY7 (40) »0Y8(10) ,0¥9(10) 9010620) ,O¥ 11010 
1), 0Y12(10) ,DV13 C10) »DV14 (10) DV15 (10) pCYAL (20) pCV12 (10) y CV13 (20) 9C 
2714010) ,CY15 (10) pCV16(10) ,CY16(103 gDY17 (10) pCV¥17(10) 
; COMMON/CACFT/FUE,A70 
INTEGER W 
GW=W/1000. 
OrS=0. 
IF CACFT.EQ. FE ) GO TO 1 
IFCACET.EQ.A70 ) GO TO 2 
1 OGW=30. 
DV1L(1)=.013 
OY1i(2)=.002 
OY11(3)=.008 
DY11(04)=-.023 
OY11 (8) =. 064 
DY11(E)=-.421 
NY11(7)=.179 
DY11(8)=-.18 
OY12(1)=.013 
OY12¢2)=.002 
0Y12(3)=.008 
0Y12(4)=-.023 
DY12(6)=.061 
DY12(€)=-0121 
OY12(7)5.179 
DY12(8)=-00 
OY13(1)=06013 
CY13(2)=.001 
DY13(3)=.008 
| DY13(4)=-.023 
DY13(&)=.062 
DY13(E)=—.121 
DY13(7)=-222 
DY13(8)=-.081 
OY14(1)=.02 


DY14(2)=.012 
OY14(3)=-.027 
DY14(4)=.037 
DY14(5)=-.005 
DY14(€)=-.090 
i OY1467)=.293 
DY14(8)=-0195 


OY15(1)=.02 
0Y15(2)=.012 
0Y15(3)=-.,027 
0Y15(4)=.037 
OY15(5)=-.005 
0Y15(€)=-.090 
0Y15(7)=.495 
0Y15(8)=-1.81 
0V16(1)=.035 
0Y16(2)=-.008 
uyY7i6(3)=-.001 
0Y16(04)=.005 


0V16(5)=.045 
OYV16 (6) =-.175 
OY16(7)=.881 
OY 16 (8) =-4.095 
0Y17(1)=.046 
OY17 (2) ==.03 
O0Y17(3)=.032 
O0V17 (4)=-.039 
OY17(5)=.1 
OV17(€)e-.141 
0V17(7)=.258 
OY17 (8)=-1.383 
CY11 (1) =.4% 
CY¥11(2)=.05 
Cvi1¢323=.05 
CV¥1104)=-.01 
CY¥11(53=-+.36 
CYi1(€)=1.2310 
CY¥1107)=-2.76 
CY¥11(8)=4.,710 
CY12¢1)0=04 
CY12(2)=.05 
C¥12(3)=.05 
CVi2(4)=-.01 
CY¥12 (5) =-.36 
CY12(CE)=1.31 
CY12(7)=-2.76 
CY¥12¢8)24.01 
CY13(1)=.4% 
CV¥13 (2)=,05 
CY12¢02)=.05 
CVi3(4)=-.01 
CY¥13(5)=-.36 
CY13(EV=RAL-31 
CY¥13(7)=-3.18 
CY¥13(8)=4,72 
CY14(1)=,.35 
CY14(2)=-.01 
CY14(4)=.24 
CY14 (5) =-41.86 
CY14 (€)=2.51 
CY14 (7) =-3.98 
CY¥14(8)=.03 
CY15(1)=.210 
CY¥15(2)=,.230 
CY15 (3) =-.43 
CY15(4)=.8 
CY15(5)=-1.76 
CY15 (€)=2,35 
CY15(7)=-8.54 
CY15 (8)=29.79 
CY16(1)=-02 
CY16(2)=.7 
CY16 (3) =-.85 


ey 
ad 


Saye enon 


CY16(4)=1.35 
CY¥16(5)=-2.75 
CY¥1€(€)=5.15 
CV16(7)2-17,9 
CY16(8)=87.65 
CY17 (1) 0h 
CY17 (2d=161 
CV17 (2) =-1.45 
CY17(04)=2.15 
CY17(5)=-3,75 
CY17(E€)=3 617 
CV¥17(7)22-36 
CYi7(@)=.21 
OY7(1)=0.0 
0Y7(2)=0.0 
OY7(3)21.6 
DY7(4)=-3,. 
0Y7¢(5)=6,.5 
OY7(6) 2-12. 
OY7(07)9=20.5 
OYe(t)=1. 
0Y8(2)=0.0 
DV8(3)=-16 
DY8(4)=4, 
OY8(5)=-10. 
OY8(6)=21. 
OY8(7)=-b1. 
OY9(1)=1. 
DY9(2)=0.0 
0Y9(3)=1.5 
OY904)=-2.5 
OYe(5)=21.5 
OY9(E) =2.5 
0Y9(7)=-6. 
OY1i0(1)=1. 
CY10(2)=2. 
OV10(3)=- 2. 
OY10(4)=2, 

\ DY10(5)=-3, 
OY10(€)=5. 
OY10(7)=-4, 
St1=0.0 
CONST1=0.0 
SL2=4. 
SL3=4, 
SL4=5, 
SL5=5, 

KK1=8 
KK2=7 
OX1=2. 
DX2=4, 
OX3=E€. 
OOT=20. 
60 TO 200 
2 DGH=20. 
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OY5(1)=.02 
0Y5¢2)=0.0 
0Y5(3)=0.0 
OY5(4)=0.0 
OY5(5)=.04 
DY5(6)=-,17 
OY5(7)=.54 
OY5(6)=-1,53 
OY6(1)=.25 
OY6¢€2)=.15 
OYE(3) =-. 30 
OY6(4)=.55 
NY6(5)=-,85 
OYE(6)=1.45 
OYE(7)=-2,5 
DYEC8)=4.25 
DY7(1)=2.0 
OY7(2)=0.0 
OV703)z-1.5 
OY7(4)=3,5 
OY7(5)=-5,75 
OY7(6)=8, 
OY8(1)=2.0 
O0Y38(2)=1.5 
BYS(3)=-1, 
NYBC4)=4,5 
OV8(5)=-5, 
OY8(6) =-4,5 
OY9(1)=2, 
0Y9(02)=3,.5 
DY9(3)25, 
CY¥904)="13,5 
OYT05)=44, 
OY9(6) =-118.5 
OV10(Cid=2, 
0Y10¢€2)=11.5 
OY10(3)=-13, 
0Y10(€4)=18,5 
OV10(5)=44, 
OY10(€)=-118.5 
CONST1=.15 
SLi=.01 
SL2=12.5 
SL3=12.5 

* SL4=12.5 
SL5=12,.5 
KKi=6 
KK2=6 
OX1=2. 
OxX2=2.5 
OX3=3, 
OOr=50. 
GO TO 13 

200 IF(GH.GT.35. ) GO TO 3 
KK1=8 
s 
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100 


102 


103 


700 


00 100 I=4,KKi 
OYS¢(TsO0YV12 07) 
OY6CT)=CYI1C7) 

GO TO 13 

TF(GH.GT.40. GO TN & 
KKi=6 

00 101 Y=1,KKi 

OYS(T) =0YV12¢7) 
OYECT)=CY12¢€T) 

GO 70 13 

IF(GW.GT.45. 9 GO TO S 
KK1=8 

00 102 T=1y)KKt 

OYS(T =0Vi3¢7) 
OY6CT=CV1I3CT) 

GO TO 13 

IF(GW.GT.50. » CO TO 6 
KKi=8 

00 102 T=1,KK1 
OYS(T)=O0Y246(T) 
OYECT=CY14(T) 

GO To 13 

IF (GW.GT.55.) GO TO 7 
KK1=8 

00 104 T=1,KK1 
OVS¢(IT)=0Y1S5(T) 

OYE CT=CYV15¢1) 

60 TO 13 

TF (GW.GT.€0.) GO TO 8 
KKi=8 

OO 105 T=1,KK1 

OYS(T) =O0Y416(T) 
OY6(T)=CY16¢T) 

GO TO 13 

KKi=8 

00 10€ T=41,KK1 
OYS(1)=0Y17¢1) 
OYECT)eCv17¢T) 

CO=1. 

U=ALT/5000. 
IFCACFT.EQeA70) UtU-i. 
CO 700 [=1,KKi 
CO=CO*U/FLOAT(T) 

U=U-1te 
SLi=DY5(1T)*CO+*SL1 . 
CONSTL=0YEC(I) *CO +CONST1 
T=SLi*°(GH-OGW) *#CONST41 
TF CACFT EQeFUE SANDS (Tol Tee eOReToGlei2e4)) RETURN A 
ITF CACFET -ECeATD AND es (Tol TeDe cORetoGTo5e75)3} RETURN A 
CO =1. 

U=DI/00! 

DO 701 T=1,KK2 
CO=U*FCC/FLOAT(T) 

UzU-i.e 

SL2=O0Y7¢(1)*CO +St2 
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RzT*St2 
IFC(T.LE.OX2) GO TO 801 
CO z1. 
V=O1/001 
00 702 T=1,kKK2 
Co=CO*U/FLOAT(T) 
UstU-1. 
702 SLI=OVBC(TI)®S04SL3 
RESL2°CXL+SLI*(T<OX1) 
IFC(T.LE.OxX2) GO TO 801 
C0=1. 
U=OI/COT 
00 703 T=i,KK2 
CO=COFU/FLOAT(T) 
U=U-1. 
703 SL&=OVOCTI*CO*SL& 
ReSLZ°OX1 *¢SL3*(OX2-0K1) *+St4*( T-0x2) 
IF (T.LE.OX39 GO TO 8601 
Cozi. 
v=01/001 
CO 708 IzigKK2 
CO=COFUs/FLOAT(T) 
UzU-1. 
708 SLS=OV10¢1)9"CC + SLS 
R=eSL2*Ox1 #SL3*(OX2-0K1) +SLGF(OX3 -0X2) #SL5%(T-0Xx3) 
801 IFCACFT.EC.FOESANDReGT.74.) RETURN A 
TFCACFTEGeFUEDAND. 2eGTo (7 .279T4+148.2)) RETURN A 
TIF CACFT EG eA7D.ANDeRoGT21000) RETURN A 
803 OseR 
GW=GW"1000. 
300 RETURN 
ENO 
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SUBROLTINE SLFCACFT SALT, Wy CI, FUEL) 
OIMENSTON DY5(10) ,076(10),0Y7(10) ,DY¥8 (10) .0¥9(10) ,0¥11(10) ,OY12¢10 
1) ,0V127(10),0¥14 (10), 0V15010) ,CV¥11 0620) »CV¥12(¢10) gCV13(10) »CY14(10) 9C 
2¥15010)50Y16(10) »CY¥16(10) ,OVL7 (10) pCV¥17 610) 
COMMCA/CACFI/F4E,A70 
INTEGER W,FUEL 
GW=W/1000. 
FUEL=Q. 
TFCACFT.EQ.FSE ) GO TO 1 
TFCACFT.EC.A70 ) GO TO 2 
1 OGw=30. 
0Y11(1)=.013 
OY11(2)=.002 : 
DY11¢2)2-.005 
OViL(4)=.021 
OY11 (6) =-.06 
OY 11 (6920137 
' OY14(7) =-.267 
0Y11(8)=.466 
0Y12¢1)=.013 
i OY12¢2)2.002 
| CY12¢3)=-.005 
0Y12(4)=.021 
0Y12(5)=-.06 
| CYL2 (€)2.137 
0Y12(7)=-.267 
0Y12(8)=.505 
0Y13(1)=.013 
0Y13(2)=.002 
0Y13(3)=.007 
DY13(4)=-.027 
| CY13(5)=.072 


DY13(E)=".154 
0Y13(7)=.287 
DY13(8)=-,3368 
DY14(1)=.020 
| CY14(2)=-.002 
| DY14(2)==,.002 
DY14(4)=.044% 
DY14(5)=-,038 
OY14(E6)=.089 
DY14(7)==,173 
DY14(8)=.388 
| DY15(4)=.02 
DY15(2)=-.002 ‘ 
0Y15(3)=.009 
DY15(4)=-,030 
OY15(5)=.094% 
CY15(E6)==5242 
DY15(7)=.601 
DY15(8)=-1.46 
CY16(1)=.022 
DY16(2)=.001 
1¥16(3)=-.006 
| 0Y16(4)=.016 
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OY 16(5)2-.026 
0Y1€6(6)=.059 
0Y16(7)=-.148 
OY16(8)=,340 
0V17(1)=.022 
OY17(02)=.001 
0Y17(3)=.005 
OY17 (4) =-,028 
OY17+5)=.1000 
OY17(6)=-.197 
0Y17(7)=.153 
OY17 (8)=.508 
CY1101)=.530 
CY11(2)=-.08 
CY11(°3)=.04 
CY11(04)=,08 
C¥11(05)2-,45 
CY11(€3=1.32 
CV11(7)=-2.94 
CY¥11(8)=5,.42 
CY12(1)=.53 
CV¥12(2)=-.08 
CV12¢3)=.04 
CY12(4)=.08 
C¥12 (5) =-.45 
C¥12(€)=1.32 
CY¥12(7)=-2.94 
CY12(6)=5.27 
CV¥13 (1) 2.530 
CY13(2)=-.08 
CV¥13(3)=-.07 
CY13(4)=.52 
CY13(E)=-1.67 
CVY13(€)=4.04 
CY13(7)=-8.13 
CY¥13(8)=12.83 
CY¥14(1)=.45 
CY14(2)=-.04 
CY14(3)=,.05 
CY14(4)=-.03 
CY¥14(5)=-,12 
CY14(E)=.,47 
CY¥14(7) 5-12.38 
CY14(8)=3.29 
CY15(1)=,45 
CY15(2)=-.04 
CY15 (3) e-.22 
CY15 (4)=1.05 
CY15 (5)=-3,15 
CY¥1566)=7045 
CY15(7)=-16.26 
CY15(8)=3€.09 
CY16(1)=.32 
CY16(2)=.01 
CY16(3)=.02 
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CY1E(4)=-,22 

CY16(5)=.95 

CY16(E)=-3.6 

CY16(7)=1049 

CY16(8)=-26,29 

CY17(41)=,532 

CY17(2)=.014 

CY17(3)=-026 

CY17(4)=.82 

CY17(5)=-2.35 

CY17(6)=3.88 

CY17(7)=,54 

CY17(8)==26.29 

CY7(1) 2045 

DY7(2)=-,25 

DY7(3)=,25 

OY7(4)=-,05 

DY7LE} 2-435 

DY7(6)=.65 

DY7(7)=.05 

DY8¢1)2.5 

DYA(2)=.25 

NY8(3)2-.5 

DY8(4)=,.75 

OY8(5)=-1, 
0Y8(6) 21475 
DY8(7) =-347 

0Y9(41)=1.25 
| DY9(2)=-4, 

DYSi3)=2. 

DY9(4) ==3,425 
| DY9(5)=4,75 

OY9(E)==3, 
| OY9(7) =-4.5 
: S$ti=0.0 
CONST1=0-0 
SL2=1.75 
SL3=1.75 
SU4=1675 
OX1i= 2.5 AT cre 4 sda rg angst 188 ARSENE AEA EIN, 
DX2= 4. ae! 
KK1=8 
KK2=7 
ODI=20. 
GO TO 200 
' 2 OGH=20. 


OY5(1)=.02 
OY5(2)=.01 
OY5¢€3)=-.02 


wane e 


0Y5(4)=.04 
0Y5(5)=-.08 
0Y5(6)=.2 
OYE(7)2-.51 
OV5(8)=1.20 
0Y6(1)=.45 
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OY6C2)=e1 
OY6(3)=-22 
OYE(4)=,2 
Ov6ét® )=-.6 
O0Y6(€)=1.55 
OYE(7)=-3.9 
OYVEC8dzSet 
OY7(1)=.16 
0Y7(2)=.08 
OY7(3)=-.18 
OY7(4)=.28 
OY7(5)=-.360 
OY7(6)=.4 
0Y68(1)=.16 
OY8(2t=.24 
OY8(3)=-.29 
OY8¢(4)=1.34 
0Y8(5)=-4.59 
OYS(EdF11.24 
OY9(1)=,4E 
OY9(2)=.04 
0Y9(3)=.36 
OY9(4) =<. 76 
0Y9(5)=2.€6 
DY9(6)=-7.56 
~ CONST1=.4 
cil=.04 
SL2+1 244 
SL3=1.44 
SL4=1.44 
OX1=2.5 
Ox2=3.5 
KK1=8 

KK2=€ 
OnT=50. 

GO TC 13 

200 TF(GW.GT 2356 


tr ee 


» GO TO 3 


aimed PD rarest OLED i OL OOOO, ad 
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KK1=8 
O00 100 T=1,kKi 
OYS(1T)20V1101T) 
100 OY6(T)zCY1101) 
GO TC 13 
3 IF(GH.GT.490. 9 GO TO 
KK1=8 
OO 101 T=1,KK2 
OY5(T) =OY120T) 
101 OVECT)=CY12¢1) 
GO TO 13 
& TF(GR.GT.&S. 1 60 TO 5 
KK1=8 
OO 102 T=i,KKi 
OYS¢(T)=O0Y13€1) 
102 DY6CTI=CYISCTS 
GO fc 13 
ITF(GW.GT.506 » GC TC 6 


wi 


KK1=8 
OO 203 T=41,kK1 
VS5SCTI=DV14CT) 
103 OYVECT)=CY14(T) 
GO TO 13 
6 IF(GH.GT.55.) GO TO 7 
* KKi=8 
00 104 Isi,kKi 
OYSCT)=O0Y15CT) 
104 OV6CT)=CY¥15(Y) 
GO TO 13 
7 TF(GW.GT.E0.) GO TO 8 
KKiz4& 
00 105 [=1,KK1 
OYS5CT) =OY26C1) 
105 DY6CT)=CY16CT) 
GO TO 13 
8 KKi=6 
DO 106 T=t,KKt 
OYS(I)=O0Y17¢0T) 
106 OY6(T)scyvi7(1) 

13 CO=1. 

U=ALT/5000. 
IFCACFT.EQ.A70) U=U-1 
00 700 I=1,KK1 
CO=CO*U/FLOAT(TY) 
U=sU-1. 
SL1=0Y5(1)*CO+SL1 

700 CONSTL=CYE(I)*CO +CONST1I 
T=SLi*(GW-OGH) +CONST1 
CO =1. 

U=0I/00TI 

CO 701 T=1,kK2 
CO=CO*L/T 
UsUe-1, 

70% SL2=OY7¢(1)*CO 4+S5L2 
R=T*SL2 : 
IF(T.LE.0X1) GO TO 8601 
CO =1. 

U=DI/001 

CO 702 T=1,KK2 
CO=U*CO/FLOAT (TY) 
UsU-1. 

702 SL3I=O0Y8CI)*CO+SL3 
R=SL2*OX1+SL3*(T-0X1) 
IFC(TSLE.CX2) GO TO 801 
CO=1. 

U=0IF/00TF 

DO 703 T=1,KK2 
CO=CO*U/FLOAT C(I) 
U=U-1. 

703 SL&=0Y9(I)*CO+SL4 

R=SL2*CX1+SL2*(OX2-DX1) +5047 7 T+OX2) 
801 CONTINUE 
FUELER*100. +.5 
j 
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GH=GW"1000. 
300 RETURN 
ENO 
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‘ SUBROUTINE CROCACFT,ALT, WyD0I, TNMPPF) 
; DIMENSION 0Y1(10),0Y2(10), OY3(10) ,D0Y4(10) 
: COMMON/CACFT/F4E,A70 
{ INTEGER W 
GW=H/1000. 

TNMPPF=0, 
YFCACFT.EQ.F4E ) GO TO 1 
IFCACFT.EQ.A70 ) GO TO 2 
41 OGW=30. 
OYiC1)-001 
DY41(2)=0.0 
0Y1(3)=0.0 
OYi(4)=-.01 
0Y1(5)=.03 
D0Y116) =-.06 
DY1(7)=.08 
0Y1(8)=-.06 
0Y2(1)=.47 
DY2(2)=,23 
OY2(3)=-.2 
OY2(4)=.09 
0Y2(5)=.2 
DY2(6) 2-272 
DY2(7)=1.27 
OY2(8) 9102 
OY3{(1)=-.002 
OY3(2)=.001 
DY3(3)=-.001 
DY3(4) =.0005 
0Y3(5)=.0013 
DY3(6)=-.0053 
OY3(7)=.0123 
O0Y4(1)2-.002 
DY4(2)=0.0 
DY4(3)=, 001 
DY4(4)==.002 
0Y4(5)=.002 
DY4(6)=.004 
OY4(7)=-.01 
SLi=-.02 
CONST1=2. 
SL2=.016 
CONST 2=.028 
KK1=8 
KK2=7 
DOI=20. 
GO TO 13 
2 OGw=15, 
0Y1(1)=0.0 
DY1(2)=-.01 
0Y1(3)=.02 
DY1(4)=-.04 
DY1(5)=.06 
1 DY1(6)=-.06 
DY1(7)=-.01 
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13 


700 


701 


DYi(8)=.22 
0Y1(9)=+.65 
OY2(1)=.386 
OY2(2)=.02 
OY2(3)=-+.902 
OY2(4)=,12 
OY2(5}=-,32 
0Y2(6)=.67 
OY2(7)=~-1,22 
OY2(8)= 1.97 
OY¥2(9)2-3.17 
OY3(1)==.01 
OY3(2)=. 0045 
OY3(3)=-+,0005 
OY3(4)2-,0056 
OY¥3(5)=.0167 
0Y3(6)==+.0357 
OY4(1)=0.0 
OY4(2)=-.01 
0Y4(3)=,02 
OY4(4)=-, 035 
0Y4(5)=.058 
OV4( 6) =-,993 
CONST1=1.12 
SL1=-.02 

CONST 2=.05 
SL2=.05 

OOT=50. 

KK1=9 

KK2=6 

CO=1. 

CO1i= 1. 
U=ALT/5000. 

00 700 I=1,KK1 
CO=CO*U/FLOAT(Y) 
CO1=CO1*U/FLOAT(I) 
SLiI=Z0V1(T)*CO4SL4 
CONST 1=DY2(1I) *COL+CONST4 
U=U-1, 

CONTINUE 

CoO=1, 

COi=1. 

U=OI/00F 

00 701 I=1,KK2 
CO=CO*U/FLOAT (TY) 
CO1=CO1*U/FLOAT (I) 
SL2=D0Y3(I)*CO +SL2 
CONST 2=DY4(1) *C04 4CONST2 
U=U-14,. 

CONT INUE 

T=(SL1* (GH-0GW) +CONST1) 
TNMPPF= SL2*(SL1*(GH-9GH) +CONST4) + CONST2 
GW=GHW*1000. 

RETURN 

ENO 
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! INPUT LIST 
COLUMN VARIABLE FORMAT 
TYPE ONE CARD: 
1-10 ACFT A Aircraft: Left-hand justify 
11-15 ALTG F Outbound altitude: feet 
16-20 ALTR F Altitude returning: feet 
21-25 ALTAT F Altitude in target area: feet 
30 Cc ! If C equals one, values of parameters are 
printed out after each segment of the mission 
for each iteration to permit checkout of the 
program. 
TYPE TWO CARDS: 
1-4 FFSAT I Fuel for start and taxi: pounds 
5:8 FFTO I Fuel for takeoff: pounds 
9-12 FFDG | Fuel for outbound descent: pounds 
17-20 FFDG | Fuel for descent returning: pounds 
21-24 FFLAT | Fuel for landing and taxi: pounds 
25-28 FTBR I Fuel to be returned: pounds 
29-32 PR { Payload returned: pounds 
33-35 DIG F Drag index outbound 
36-38 DIR F Drag index returning 
41-45 FAS I Fuel at engine start: pounds 
46-50 P | Payload (any item capable of being dropped, 
ejected or jettisoned): pounds 
51-56 BAW ! Basic aircraft weight (gross weight at engine 
Start less fuel and payload): pourds 
57-60 DFSAT F Distance for start and taxi: NM 
33 
i 
| 
h 


COLUMN VARIABLE FORMAT 


61-64 DFTO LS 
65-68 DFDG F 
69-72 DFL F 
73-76 DFDR F 
77-80 DFLAT F 


Distance for takeoff: NM 

Distance foi outbound descent: NM 
Distance during loiter in target area: NM 
Distance for descent returning: NM 


Distance for landing and taxi: NM 


TYPE THREE CARD - An EOR (7/8/9) card follows the last Type Two card of a set of data 


Any number of sets of data may be used. 


APPENDIX III 


PROGRAM LIMITATIONS 
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PROGRAM LIMITATIONS 

F-4E and A-7D Aircraft 
All climbs are at military thrust. 
All cruises are at optimum Mach numbers. 
Outbound and returning altitudes: 5,000 to 40,000 feet 
Altitude in target area: 0 or 5,000 to 40,000 feet 
F-4E restrictions: 

Gross weight: 30,000 to 58,000 pounds 

Drag indicies: O to 140 
A-7D restrictions: 

Gross weight: 20,000 to 42,000 pounds 


Drag indicies: 0 to 300 
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F-4E DRIVER FOR RADIUS OF ACTION 
SUBROUTINES LISTING 
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A-7D DRIVER FOR RADIUS OF ACTION 
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TABULAR DATA FOR THE F-4E AIRCRAFT 
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43 aestTwmacr 
: A computer program using an iterative technique was developed for determining the radius of action 
of F-4E and A-70D aircraft with any configuration of external stores. The program uses perforrnance 
data from the aircraft flight manuals to calculate fuel and distance required to achieve military power 
climbs and optimum cruise fuel consumption. Required inputs to the pregram include the initial 
amount of fuel onboard and aircraft gross weight, the cruise and loiter altitudes, and the outbound 
and returning drag indictes due to aerodynamic drag of the external stores. Optional inputs that 
provide increased accuracy include fuel for angine start and taxi, and fuel and distance for takeoff 
and descent. The program calculates sequential increments of the outbound cruise portior of the 
mission profile until the fuel reserve desired at the end of mission (return to home base) is obtained; 
thus the radius of action is determined. 
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